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We report on a study of the electronic and magnetic properties of the triangular antiferromagnetic 
{Cits} single- molecule magnet, based on spin density functional theory. Our calculations show that 
the low-energy magnetic properties are correctly described by an effective three-site spin s = 1/2 
Heisenberg model, with an antiferromagnetic exchange coupling J « 5 meV. The ground state 
manifold of the model is composed of two degenerate spin 5* = 1/2 doublets of opposite chirality. 
Due to lack of inversion symmetry in the molecule these two states are coupled by an external 
electric field, even when spin-orbit interaction is absent. The spin-electric coupling can be viewed 
as originating from a modified exchange constant 5J induced by the electric field. We find that 
the calculated transition rate between the chiral states yields an effective electric dipole moment 
d = 3.38 X lO'^^C m ^ elO"*a, where a is the Cu separation. For external electric fields e ~ 10 
V/m this value corresponds to a Rabi time r ~ 1 ns and to a 5J of the order of a few jieV. 



I. INTRODUCTION 

Single-molecule magnets (SMMs) have been in- 
tensively studied in the last two decades (for a re- 
view see Ref. [l]). At low temperature these remark- 
able molecules behave in part like bulk magnets 
thanks to their very long magnetization relaxation 
time. At the same time SMMs are genuine quan- 
tum systems. They display a variety of non-trivial 
quantum effects such as the quantum tunneling of 
the magnetization)2i^ Berry phase interference?^ and 
quantum spin coherence^-. Due to their double na- 
ture, SMMs are ideal systems to investigate decoher- 
ence and the interplay between classical and quan- 
tum behavior 

From the point of view of applications, interest in 
SMMs has been in part spurred by the possibility 
that these structures could represent the ultimate 
molecular-scale limit for magnetic units in high- 
density magnetic storage materials. More recently 
SMMs have been recognized as promising building 
blocks in molecular spintronics, the emerging field 
combining spintronics and molecular electronics i^^— 
In particular, thanks to their long spin coherence 
time?^ SMMs are good candidates to realize spin- 
tronic devices that maintain, control and exploit 
quantum coherence of individual spin states. These 
devices could find important applications in the field 
of quantum information processing] ^^d^ 

One key issue in using SMMs in molecular spin- 
tronics and quantum information processing is the 
ability of switching efficiently between their different 
magnetic states. The conventional way of manip- 
ulating magnetic states is by applying an external 
magnetic field. However, this approach has signif- 



icant drawbacks when it comes to controlling mag- 
netic states at the molecular level. Quantum manip- 
ulation of SMM requires application of an external 
field at a very small spatial and temporal scale. It is, 
however, very difficult to achieve such a small scale 
manipulation using standard electron- spin control 
techniques such as electron spin resonance (ESR) 
driven by ac magnetic field. ^ 

One promising alternative to achieve control of 
magnetic states at the molecular level is to use an 
electric field instead. Typically, by using STM tips 
for example, it is possible to apply strong time- 
dependent electric fields in sub-nano regions, with 
time scales of 1 nsJ^ii^ Clearly, since electric fields 
do not couple directly to spins, it is essential to find 
efficient mechanisms for spin-electric coupling as well 
as real SMMs where this mechanism can be at play. 
In principle an electric field can interact with spins 
indirectly via the spin-orbit interaction. However, 
since the strength of the coupling scales like the vol- 
ume of the system, this mechanism is not the most 
efficient one for manipulating SMMs. 

Recently, it has been proposedi^ that in some 
molecular antiferromagnets lacking inversion sym- 
metry, such as the triangular antiferromagnetic 
{Cwa} and other odd-spin rings, an electric field 
can efficiently couple spin states through a com- 
bination of exchange and chiralilty of the spin- 
manifold ground state^-. The {Cus} molecule, while 
largaiiii^, reduces to a simple model composed of 
three identical spin s = 1/2 Cu cations coupled 
by an antiferromagnetic (Heisenberg) exchange in- 
tegration. Its ground state consists of two total-spin 
S — 1/2 doublets of opposite spin chirality, degener- 
ate in the absence of spin-orbit interaction. Accord- 



ing to an analysis based on group theory ; -^^i-^^ due 
to the lack of inversion symmetry, an electric field 
can couple states of opposite chirality through the 
dipole operator, even when spin-orbit interaction is 
absent. In the presence of an additional small dc 
magnetic field that mixes the spin states, this spin- 
electric coupling will then generate efficient electric 
transitions from one spin state to another. 

An intuitive picture of this coupling is the fol- 
lowing. Since a spin S — 1/2 triangular antiferro- 
magnet is frustrated, there exist three energetically 
degenerate antiferromagnetic spin configurations for 
= 1/2 and three for = -1/2. Both ground 
state chiral eigenstates, with a given value of Sz, are 
appropriate, equally weighted, linear combinations 
of these three frustrated spin configurations. Each of 
these three configurations, if prepared, would have a 
dipole moment with the same magnitude that points 
from the antiparallel sites to the midpoint between 
the two parallel sites. While the net dipole moment 
of the two chiral eigenstates is zero, the dipole tran- 
sition matrix element between them is not and it is 
simply related to the magnitude of the permanent 
dipole moment of the energetically degenerate frus- 
trated configurations. 

In practice the relevance of this spin-electric mech- 
anism depends on the coupling strength of the chiral 
states by the electric field, i.e. on the value of the 
dipole moment of the frustrated spin configurations. 
Theoretically this is an issue that only a micro- 
scopic calculation for the specific molecule can ad- 
dress. The main objective of this work is to calculate 
the strength of this coupling for the {Cu^} molecule 
using ab-initio methods. Our approach is based 
on Spin Density Functional Theory (SDFT), imple- 
mented in the NRLMOL codes, which has been very 
successful in describing the electronics and magnetic 
properties of Mni2-acetate and other SMMs.— 
Recently SDFT implemented in NRLMOL has been 
used in a first-principle study of quantum transport 
in a Mni2 single-electron transistor;2i 

Our results show that indeed the crucial electric- 
dipole moment is not negligible in {Cua} and it 
would correspond to characteristic Rabi times of 1 
ns in the presence of typical electric fields generated 
by STM tips. As originally suggested in Ref. [H, 
the spin-electric coupling can be interpreted as due 
to a modified exchange interaction brought about 
by the electric field. Although here we only address 
the specific case of {Cua}, our paper introduces a 
methodology that can be followed in a systematic 
study of other SMMs without inversion symmetry. 

The paper is organized as follows. In section |lT] 
we discuss the electronic and magnetic properties of 
triangular {Cua} molecule based on ab-initio calcu- 
lations and show that the low-energy quantum prop- 



erties of the molecule can be described by an effectiv? 
three-spin s = 1/2 Heisenberg model with antiferro- 
magnetic coupling. In section III Bl we review the 
underlying mechanism of spin-electric coupling in 
{Cwa} antiferromagnet, based on the effective spin 
Hamiltonian. The first-principle computation of the 
spin-electric coupling and electric dipole moment of 
{Cita} is presented in section IIVI In section IIVBI 
we discuss the effect of the electric field on the ex- 
change coupling. Finally we present the summary of 
our work in section IVl 



II. ELECTRONIC AND MAGNETIC 
PROPERTIES OF {Cu^} 

A. Microscopic description of the molecule 

The {Cwa} molecule that we are interested in 
has chemical composition Nai2[Cu3(AsW9 33)2- 
3II2O] -321120^. This molecule has been studied 
experimentally by different groupsJ^ii^ The three 
Cu^^ cations form an equilateral triangle and, as we 
show below, are the sites of three identical s = 1/2 
quantum spins. The frontier electrons on each of 
these sites have primarily d character. The bridging 
atoms consist of predominantly paired electrons and 
are only polarized to the degree that the same-spin 
states hybridize with the unpaired d-electrons on the 
Cu sites. Due to the localized nature of transition- 
metal 3d states, direct exchange stabilization due 
to parallel neighboring states is expected to be ex- 
ponentially small. Therefore, unless the frontier d- 
electrons are spatially orthogonal by symmetry to 
the d-electrons on other sites, antiferromagnetic or- 
dering between electrons on a pair of neighboring 
Cu atoms is energetically preferred due to the in- 
crease in the system's kinetic energy, induced by or- 
thogonality constraints, when neighboring states are 
parallel. 

Although the spin model of three exchange- 
coupled spin 1/2 is quite useful to understand the 
magnetic properties of the {Cwa} SMM, all the other 
atoms in the molecule are essential for its geometri- 
cal stability and for the resulting superexchange in- 
teraction among the spins at the Cu sites. A proper 
ab-initio description of the molecule must therefore 
include to a certain extent all these atoms. 

Building a suitable model of the molecule is a con- 
siderable challenge since the model molecule should 
preserve the essential physics. We have constructed 
the molecule by preserving the Da/j symmetry of the 
polyanionic part of the molecule as observed in the 
experiment jiiii^ Three of the twelve Na atoms of 
the molecule are placed at the belt region of the 
molecule. These three are the most important of all 



the Na atoms for the stability of the belt region of 
the molecule. There is some uncertainty in the posi- 
tion of the Na atoms but we have placed eight of the 
remaining nine Na atoms in a way to preserve the 
D3/1 symmetry. The last Na atom is replaced by a H 
atom and is placed at the center of the molecule to 
maintain the charge neutrality of the valance elec- 
trons. The model of the molecule used in this calcu- 
lation is shown in Fig. [TJ 




FIG. 1: Model of the {Cua} molecule with chemical 
composition NaiiH[Cu3(AsW9 33)2- 3H2O] used in this 
work. Xcrysden visualization tool^^ is used for this fig- 



We have relaxed the geometry using the a&-initio 
package NRLMOL26.27 ^^^^ 

uses a Gaussian basis 
set to solve the Kohn-Sham equations using PBE- 
GGA approximation. 2* All-electron calculations are 
performed for all elements of the molecule except 
for tungsten, for which we have used pseudo poten- 
tials. The relaxation is first performed by setting the 
net total spin of the molecule to S = 3/2 and then 
by changing the net spin to 1/2. Self-consistency is 
reached when the total energy is converged to 10^^ 
Hartree or less. 

The density of states of the molecule is shown in 
Fig. m The HOMO-LUMO gap for the majority 
spin is calculated to be about 0.78 eV and that for 
minority spin is about 0.58 eV. Although in our cal- 
culations we have used an equilateral arrangement of 
the three Cu atoms, it is found experimentally that 
the {Cus} molecule in the ground state is slightly 
distorted into an isosceles triangle. Since the cal- 
culated HOMO-LUMO gap for the equilateral con- 
figuration is relatively large, the distortion is likely 
to be due to magnetic exchange rather than to the 
Jahn- Teller effect. 

One important result of our calculations, after the 
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FIG. 2: Density of states of {Cus} molecule. HOMO- 
LUMO gaps for majority and minority spins are shown 
in the inset. 



geometry relaxation have been implemented, is that 
the ground state of the system is anti-ferromagnetic, 
with a net total spin of 1/2 in accordance with 
experiment 1^ The ground state energy is lower by 
about 5.6 meV relative to spin S = 3/2 configura- 
tion. This allows us to assign an exchange constant 
J « 5 meV to the three-site Heisenberg spin model 
mentioned above (see also next Section). 

The calculated magnetization density of the re- 
laxed molecule shows the presence of three electron- 
spin magnetic moments fii ~ 0.55/xb, i = 1, 2, 3, es- 
sentially localized at the three Cu atom sites. Note 
that the orbital moments are quenched. These re- 
sults confirm that the low-energy properties of the 
{Cus} molecule can be approximately described by 
an effective spin Hamiltonian of three spins s = 1 /2 
localized at the Cu sites. 

The exchange coupling between two Cu atoms is 
indirect and follows a superexchange pathi^ along 
Cu-O-W-O-W-O-Cu as shown in Fig.H- see the ye- 
llow line connecting the atoms. To understand this 
coupling mechanism we focus on one of the three 
CuOs complexes of the molecule (shown inside the 
circle in Fig. Because of the square-pyramidal 
Civ point-group symmetry of this complex, the dxy , 
dxz, dyz states of Cu have lower energies compared 
to the d^2_y2 and states. Moreover, our calcu- 
lation shows that the axial Cu-0 distance (2.35 A) 
in each unit is larger than the four equatorial Cu- 
O distances (1.93 A). Thus the energy of d^^ state 
is lower than d^2_y2 state and the unpaired d elec- 
tron of the Cu^+ ion resides in dx2-y2 state that is 
directed along the equatorial Cu-0 vectors. There- 
fore, the exchange coupling between two Cu atoms 
involves three O atoms and two W atoms. 

The magnetic moment calculations of the atoms 
of {Cits} molecule also support the superexchange 
path. The magnetic moments at the O and W atoms 
on this path is much smaller than at the Cu sites, 



FIG. 3: (Color online) Superexchange coupling between 
two Cu atoms. The yellow line connecting two Cu atoms 
through three O and two W atoms shows the path along 
which spin coupling between Cu atoms is mediated. The 
numbers near the atoms are the magnetic moment (in 
units of /is) of the atoms along the exchange path. 



but still 2 order of magnitude larger than at atoms 
not belonging to this path. 



B. Effective spin Hamiltonian description 

Based on the results of the a&-initio calculations, 
the low-energy properties of the {Cu^} molecule can 
be described by the following quantum spin Hamil- 
tonian 



, (1) 



1=1 



where J is the exchange parameter, D is the 
Dzyaloshinski vector and are three spins-1/2, 
located at the Cu sites. The first term in the 
Hamiltonian is an isotropic Heisenberg model. The 
geometry-relaxation and electronic-structure calcu- 
lations showed that the Cu atoms form an equilat- 
eral triangle with a very small intrinsic deformation. 
Since the atomic environment around each of the 
the three Cu-Cu bonds is the same, we take the 
three exchange constants Jis+i to be the same value, 
J. On the basis of the splitting between the ferro- 
magnetic and antiferromagnetic configurations dis- 
cussed in the previous section, J is positive and ~ 5 
meV. The second term in Eq. ([l]) is the anisotropic 
Dzyaloshinski-Moriya exchange interaction origina- 
ting from spin-orbit interaction. Its strength |Di | 
is at least one-order of magnitude smaller than the 
isotropic exchange constant J, and we will disregard 
it for the moment. 

The ground state of Eq. ^ is total spin S" = 1/2 
manifold, which can be constructed in terms of six 
degenerate spin configurations, three associated with 




FIG. 4: a) The three spin configurations of the molecule 
associated with total spin projection = +1/2. b) 
The two chiral states formed from a) with chirality -|-1 
and -1, respectively. 



Sz = +1/2 and the other three associated with 
Sz = —1/2- Fig. |4] shows the three possible spin 
configurations associated with Sz = +1/2. 

The total-spin S — 3/2 four-dimensional subspace 
has an energy of order J above the ground-state 
manifold. 

Within the S* = 1/2 ground-state manifold, we can 
construct two degenerate, linearly independent dou- 
blets. Specifically the two Sz = +1/2 states (shown 
in Fig. mb)) are 
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where uj — e~ . The quantum numbers i?+ and 
E- specify the so called handness or chirality of the 
states \E± , M) , which are eigenstates of the chirality 
operator 



Cz = -^Si • S2 X S3 , 



(4) 



with eigenvalues ±1 respectively. It is useful to in- 
troduce also the other two components of the chiral 
vector operator 



Cx 



--(Si • S2 - 2S2 • S3 + S3 • Si) , (5) 



= -82 - S3 • Si) , (6) 

and the ladder operators C± = Cx i iCy. Note 
that [Ci,Cjn] = «2e/„i„C„ and [Ci,Sjn] = 0. Here 
£imn is the Levi-Civita symbol. The ladder opera- 
tors reverse the chirality of the states: C±\Ezp, M) = 
\E±,M). Thus C behaves exactly like the operator 
S (for S — 1/2) in chiral space. 

In the microscopic description of the molecule im- 
plemented within density functional theory via the 
NRLMOL code, the chiral states defined in Eqs. ^ 
and ([3]) have to be understood as being composed 
both of a spin and an orbital part. 

We conclude this section with an observation of 
the DM interaction. As shown in Ref. llQ, the DM in- 
teraction within the 5 = 1/2 ground state manifold 
takes the simple form Hdm = ^soCzSz, where Aso 
is the effective spin orbit coupling constant. Thus 
equal-spin states of opposite chirality are split by 
2Aso. 



III. SPIN-ELECTRIC EFFECT IN {Cus} 

A. Absence of inversion symmetry and 
coupling of ground-state chiral states 



turns out that in triangular spin-1/2 antiferromag^ 
nets the coupling is non-zero. On the other hand 
in pentagon spin 1/2 antiferromagnets, the coupling 
vanishes, unless spin-orbit interaction is included 

We focus now on the spin-electric coupling of chi- 
ral states in {Cus}. In the presence of an external 
electric field e. The Hamiltonian acquires the addi- 
tional electric-dipole term Hg — J^i er^ • £ = eR • e, 
where e is the electron charge and r.^ is the coordi- 
nate of the ith electron. 

In the subspace of spin projection 5*2 = 1/2 of the 
ground-state manifold, which is invariant for the 
perturbed Hamiltonian Hq + can be expressed in 
the basis of the chiral states as 

H = Hq + H, 

{E+, +^\Ho\E+,+^) {E+, +i|H,|£;_, +i) 

{E^,+l\Ha\E^,+^) 

(7) 

A similar expression holds for the Sz ~ ^1/2 sub- 
space. The eigenvalues of H are 

£;±(£) = S±(0)±|d.e| , (8) 

2 2 

with Ef{Q) = {E±,+l\Ho\E±,+\), and the corre- 
sponding eigenstates 



The triangular spin-1/2 antiferromagnet {Cua} 
belongs to the class of antiferromagnetic rings with 
an odd number of half-integer spins. In these sys- 
tems, the lack of inversion symmetry of the molecule 
as a whole implies that the ground-state is a four- 
dimensional manifold, whose basis states \E±^Sz = 
±1/2) are characterized by the spin projection Sz = 
±1/2 and by the chirality Cz = ±1 (which we also 
label as E± ). In contrast, antiferromagnetic rings 
with an even number of spins have non-degenerate 
S = singlet ground state. According to the origi- 
nal proposal in Ref. [l6lfl9l in odd-spin rings the two 
states of opposite chirality \E±,Sz = M) can be 
coupled linearly by an external electric field, even 
in the absence of spin-orbit interaction. In order 
for electric coupling to be non-zero, other criteria 
must be satisfied. First of all, permanent electric 
dipoles dij must be present on the bridges that me- 
diate the coupling of spin and Sj. A necessary 
(although not sufficient) condition for this is that 
the superexchange bridge that magnetically couples 
Si and Sj lacks a center of inversion symmetry. Even 
when local dipole moments are present on individ- 
ual bridges, the resulting final spin-electric coupling 
between chiral states depends in a nontrivial way 
on the overall symmetry of the molecule. The best 
way to settle this issue is to carry out a system- 
atic symmetry analysis based on group theory. It 
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(9) 

Here we have introduced the electric dipole matrix 
element d, which couples states of opposite chirality 
(but with the same spin projection) 

d =(£;+, +i|eR|i?_,+i). (10) 

For the specific example of {Cu^} molecule only 
the matrix elements of X and Y components of R 
are nonzero and 

{E+,+^\eX\E^,+^)=z{E+,+^\eY\E^,+^) = -| , 

(11) 

where d= |d|. 

The matrix element in Eq. (llOp is the key quantity 
in the spin-electric coupling mechanism. Substitut- 
ing the expressions for the chiral states from Eqs. ([2]) 
and using the orthogonality of spin states we obtain 

d = i((*;nleR|*;n> +^(*ntleR|*nT> + 

UJ^^ni\eIi\^^^i,)) . (12) 

Evaluating the dipole matrix element between two 
states of opposite chirality is therefore equivalent to 



calculating the dipole moment of each of the three 
spin configurations. This matrix element determines 
the strength of spin-electric coupling and we are pri- 
marily interested in calculating this quantity by a6- 
initio methods. 

Finally, note that all the matrix elements of 
the electric dipole operator eR are identically 
zero in the S = 3/2 subspace. This is obvi- 
ous since (4'^^^|eR|4'^^^) and i((vI'^^^|eR|vI'^^^) 4- 
(^'^^^|eR|^'^^^) + (^'^^;|eR|^'^^;)) are both zero by 
symmetry. We will confirm this result by direct a6- 
initio calculations. 



B. Effective spin Hamiltonian description 

The effect of the electric field on the the low- 
energy spectrum of {Cus} can be recast in the form 
of the effective spin model introduced in Sec. IIIBI 
Since the electric dipole operator has nonzero matrix 
elements only in the ground-state manifold, where 
it couples states with equal spin components and 
opposite chirality, we expect that the spin-electric 
Hamiltonian iJ^ can be rewritten as a linear com- 
bination of the ladder operators C± . By comparing 
the matrix elements of given in Eq. ([TU| and (fTTT) 
with the action of C± on the chiral states, one can 
show thatii 

iff = A£'.c,|, (13) 

where e' = R^{(j)){7TT/6 - 20)e, with i?((/)) being the 
matrix representing a rotation by an angle 4> around 
the z-axis, and 9 being the angle between the in- 
plane component ey of the electric field and the bond 
Si — S2. By using Eq. (O and ^ we can now rewrite 
C|| — {C'x,Cy) in term of spin-operators and we 
obtaii>i^ 

3 

iff = ^(5J,,+i(e)s, -sj , (14) 

i 

where the modified exchange parameters take the 
fornii^ 

Ad 27r 
SJu+iie)^^\e\\\cosi—i + e). (15) 

This expression of the effective electric-dipole 
Hamiltonian suggests a transparent physical in- 
terpretation of the spin-electric couping mecha- 
nism.— An external electric field changes the 
charge distribution of the {Cms} molecule which, 
in turn, changes the exchange interaction between 
neighboring atoms. Since the modified exchange in- 
teraction does not commute with Hq, it can cause 



transitions between chiral states within the grount^ 
state manifold. 

In Eq. ([15]), e|| is the projection of electric field 
on the Cm3 plane (in our case ey — e), i=l and 
9 = 30° is the angle between e and the line joining 
Cui and Cu2- Finally, note that Eqs. (fT4)) and ((T5|) 
provide an estimate of the dependence of the ground- 
state energy as function of the electric field. Since 
in the absence of spin-orbit coupling the electric- 
dipole Hamiltonian has zero matrix elements in the 
5 = 3/2 subspace, Eq. P3)) gives us an estimate of 
the dependence of the exchange constant J (propor- 
tional to the splitting between the S' = 1/2 ground 
state and S = 3/2 excited state) on e. 



IV. ^B-INITIO EVALUATION OF THE 
SPIN-ELECTRIC COUPLING 

A. Calculation of the electric dipole moment 

To construct the chiral states of the full {Cu^} 
molecule, we have calculated the ground state of the 
molecule for different spin configurations, as shown 
in Fig. m Although there are two doublets of chiral 
states for the triangular arrangement of three spin 
1/2 atoms, in this calculation we have used only one 
doublet associated with the spin projection -1-1/2, 
since we are interested in coupling between states of 
opposite chirality with the same spin projection. 

To study the spin-electric effect we have applied 
an external field along the perpendicular bisectors 
between positions 2 and 3 of the Cu^ triangle shown 
in Fig. [5l and have calculated the corresponding 
ground state energy self-consistently for different 
spin configurations. We have kept the direction of 
the field relative to coordinate axes fixed, and have 
changed the orientation of the spins at the Cu atoms 
to generate the three possible spin configurations of 
the {Cua} molecule. 




FIG. 5: The direction of the applied electric fields used 
in this calculation. 



Our calculations show that {Cu^} molecule in 
the spin Sz = 3/2 state does not have any per- 
manent electric-dipole moment. On the other hand 
each of the three frustrated spin Sz = 1/2 con- 
figurations have a small permanent (i.e, zero-field) 
dipole moment, as expected from the general dis- 
cussion of Sec. III. The three moments have all the 
same magnitude but their directions are along the 
perpendicular bisector of Cu^ triangle and between 
two Cu atoms with parallel spin alignments. The 
relative orientations of these moments along with 
components are shown in Fig. [51 The fact that the 
Sz =3/2 state does not have permanent dipole mo- 
ment whereas Sz = 1/2 states do, suggests that the 
dipole moments are solely due to spin effects. 
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FIG. 6: Dipole moments of three spin configurations 
and their relative angles. Pitt ~ (^4-ttk-f^l''''4.tt)' 
Pnt = (^ntleR-l^nt) and pftj, = (*tn|eR|*ft-l.) are 
the moments corresponding to the spin configurations of 
Fig.lli. 

In the presence of an electric field the energies 
of the {Cms} molecule are slightly lower when field 
is between two Cu atoms with parallel spins than 
for the other two spin configurations, where the 
field is between two Cu atoms with anti-parallel 
spin alignments. This difference in energy is due to 
the direction of permanent moment relative to the 
induced moment. We have calculated the perma- 
nent dipole-moment of the ground state spin config- 
uration by fitting the dependence of energy of one 
of the Sz ~ 1/2 spin configurations with external 
field, as shown in Fig. [T] The calculated values 
of the permanent dipole moment and polarizability 
of {Cus} molecule are p = 4.77 x 10~^^C.m and 
a = 1.025 X IQ-^^Cm'^/V, respectively. Although 
there is no experimental value of polarizability avail- 
able for {Cua}, polarizabilities within DFT calcula- 
tions are generally accurate to 1-3 percent. 

The value p extracted from this fitting is consis- 



tent with the direct calculation of the electric dipol? 
moment of the three spin configurations at zero field, 
implemented in the NRLMOL. 
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FIG. 7: Electric field dependence of the energy for one 
of the three spin Sz = 1/2 spin configurations. The plot 
for the other two configurations is very similar and the 
fitting yields essentially the same values of p and a. 

To calculate the matrix element d given in 
Eq. 1121), we substitute the components of the mo- 
ments for the different spin configurations of Fig. [51 

d = ip[(l + oj cos f3 + oj'^ cos 2^)x 
+ {uj sin/3 + a;^ sin2^)y] 

= |(x + zy). (16) 

The magnitude of the dipole coupling in {Cu^} 
molecule is, therefore 

d=-^= 3.38 X IQ-^^C.m . (17) 
v2 

The efficiency of the {Cu^} molecule as a switch- 
ing device depends on how fast an electric field can 
generate transitions from one chiral state to the 
other. The characteristic (Rabi) time for transitions 
between the two chiral states is given by 



Here, h is Planck constant, d is the dipole matrix 
element between states of different chirality given by 
Eq. pTj) . and e is the external electric field. Fig. HI 
shows the dependence of the Rabi time on external 
field, with the maximum value of ~ 50 ns for a field 
£ ~ 5 X 10^ V/m. For larger fields of the order 
of ~ 10^ V/m, easily attainable in the vicinity of a 
STM tip, the Rabi time is of the order of 1 ns, which 
is considered to be a relatively fast control-time in 
quantum information processing. 
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FIG. 8: Electric field dependence of the Rabi time for 
quantum transitions between the two (ground-state. 



B. Modification of the exchange coupling in an 
electric field 

To calculate the dependence of the exchange cou- 
pling J on the electric field, we need to determine 
how the spin S = 1/2 ground state and the spin 
5 = 3/2 excited state depend on the field. We de- 
fine the exchange energy J{e) as the difference 



(19) 



where Ei (e) and E^ (e) are the energies of the S = 
1/2 ground-state and of the spin S = 3/2 excited 
state respectively in the presence of an electric field. 

Based on on our discussion of Sec. IIIII A [see 
Eq. ([5])], the energy of the 5=1/2 chiral ground- 
state manifold and the 5 = 3/2 excited state vs e 
are shown schematically in Fig.|9l where we have dis- 
regarded the quadratic dependence of both Es and 
El on the field due to the induced electric dipole 
moment. 




FIG. 9: Schematic electric-field dependence of the ener- 
gies of the 5 = 1/2 chiral states and spin S = 3/2 excited 
state, and the exchange energy J defined in Eq. (|19|l . 

The calculation of the electric-field-modified ex- 
change parameter using first-principle methods is 



not completely straightforward, since the SDFT ca^ 
culations done within NRLMOL allow us to calcu- 
late the energy of a given spin configuration, whereas 
the (chiral) ground-state is a linear combination of 
three possible spin configurations. However, we can 
get an estimate of the dependence of J on e by ap- 
proximating 



E,{e) w ai(^';^tl-H"cFT(£)|*in) 

2 

+al{-9ut\HDFTie)\-^nt) 
+al{'^ni\HDFTie)\-^ni) 



The coefficients a's can be obtained by expanding 
X7(£)^ in Eq. ([9]) in terms of the spin configura- 
tions, which leads to 



[(1 - r)^'^^-t- + - w^r)*-^^^ 



1 

+{uj^ - wr)^'^^;] 



where r=-^^if'- = ^ *)> ^'•^'^ given choice of 

the electric field direction. 
Therefore, 




The energies E^^^ and E^^^ are the same because 
of symmetry. Since the difference between E^^^ and 
E^^^ is very small and near the accuracy limit of our 
calculations, we further approximate E^^-^ w ^tit- 

The exchange parameter J becomes 



J(e) w E^^^{e) - Ei^^{e) , 



(22) 



with i?-f^-f(e) = £'3 ((e). 

In Fig. [10] we plot the electric-field-induced vari- 
ation of the exchange energy 6J{e) = J{e) — J(0) 
vs. e. The result for 6J obtained by evaluating 
Eq. ([22]) with SDFT is shown by the red curve. For 
this part of the calculations the convergence crite- 
rion has been increased up to 10~* Hartree. We 
can see that the dependence of J on electric field 
is quite small, and 5 J is in the fieV range for elec- 
tric fields e = (1 - 10) X 10^ V/m. These energies 
are not far form the accuracy limit of our numerical 
calculations, which is the reason of the fluctuations 
seen in the plot. Nevertheless the overall trend is 



an increase of SJ(e) with e, which is approximately 
Unear at low fields. Note that the SDFT evaluations 
of E^-]-^{e) and E^^^{s) contain a quadratic contribu- 
tion in e but this nearly cancels at small fields when 
computing SJ, and it becomes appreciable only at 
£ > 5 X 10"^ V/m. 

The blue line in Fig. [TU] shows the dependence of 
dJ on e given by the prefactor of the cosine func- 
tion in Eq. (jlSp , which was derived within the spin 
Hamiltonian formalism. When plotting Eq. (jl5p we 
have used the value of d extracted from our first- 
principle calculations. Comparing the two curves, 
we note that, apart form the fluctuations in the nu- 
merical result mentioned above, the theoretical and 
numerical values for SJ are consistent, and both pro- 
cedures predict an overall increase of 6 J with electric 
field. 




Electric field, E (V/m) ^^^7 

FIG. 10: Electric-field dependence of the variation of 
the exchange energy SJ{e) = J{e) — J(0) induced by the 
field. The red curve is the first-principle result obtained 
by evaluating Eq. (I22p and the dashed black curve is 
the quadratic fit of 5J{e). The blue curve is a plot of 
Eq. I|15p with the numerical value of d extracted from 
the first-principle calculations. 



V. SUMMARY 

In this paper we have carried out a first- 
principle study of the spin-electric coupling in single- 
molecule magnets (SMMs) without inversion sym- 
metry. Specifically, we have analyzed the clear-cut 
case of the {Cu^} triangular antiferromagnet where, 
because of spin frustration, the ground-state con- 
sists of two generate spin 1/2 doublets of opposite 
chirality. Theory predict o^^'^^ that an electric field 
can couple these states, even when spin-orbit inter- 
action is absent. The main goal of our work has been 
to compute how strong this coupling is. 

Our calculations of the electronic structure of the 
{Cua} molecule show that the spin magnetic mo- 



ments are localized at the three Cu atom sites of the 
molecule. The magnetic properties of the molecule 
are correctly described by a triangular spin s = 1/2 
Heisenberg antiferromagnet, with an exchange cou- 
pling J of the order of 5 meV that separates the 
energies of the spin-S* = 1/2 ground-state many- 
fold and the spin-S* — 3/2 excited states. In agree- 
ment with theoretical predictions fi^i^ we find that 
an electric field couples the two ground-state dou- 
blets of opposite chirality, even when spin-orbit in- 
teraction is absent. The strength of the coupling 
is linear in the field and proportional to the perma- 
nent electric dipole moment d of the three frustrated 
spin configurations. The calculations yield a value 
of d « 4x 10~^^C m « elO~'^a for {Cus}, where a is 
the Cu atom separation. Corresponding Rabi times 
for electric-field-induced transitions between chiral 
states can be as short as 1 ns, for electric fields of 
the order of 10* V/m, which are easily produced 
by a nearby STM tip. Thus this spin-electric cou- 
pling mechanism is of potential interest for the use 
of single-molecule magnets in quantum information 
processing as fast switching devices. 

Our calculations also indicate that the presence 
of an external electric field modifies the exchange 
constant J. Typically the electric field increases J, 
although the energy scale of this change is in the 
^eV range for typical STM-generated electric fields. 
Thus for this specific antiferromagnetic SMM, the 
electric field cannot trigger directly a level crossing 
between magnetic states with different total spin, as 
suggested recently for other SMMsi ^"'^^ 

This work shows that a microscopic investiga- 
tion of the spin-electric coupling using the NRL- 
MOL first-principle code is feasible, and can sys- 
tematically implemented for a large class of SMMs 
which lack inversion symmetry. In this paper we 
have disregarded the effect of spin-orbit interaction 
and external magnetic field. The spin-orbit inter- 
action strength is small compared to the exchange 
coupling J. In the case of {Cu^} it simply intro- 
duces a small splitting between the chiral states, but 
is not expected to infiuence significantly the spin- 
electric coupling. However in other antiferromag- 
netic rings with an odd number of spins spin-orbit 
interaction is essential for the very existence of the 
coupling mechanismii^. Work to include both spin- 
orbit interactions and an external magnetic field is 
in progress. Together with the group-theory analysis 
presented in Ref. [19|, these studies will be a consider- 
able help in guiding future experiments and selecting 
the most promising SMMs for applications in quan- 
tum information processing and nanospintronics. 
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